Understanding correlated electron systems by 
a classification of Mott insulators 



Subir Sachdev 

Department of Physics, Yale University, 
P.O. Box 208120, New Haven, CT 06520-8120, USA 



Abstract 

This article surveys the physics of systems proximate to Mott insulators, and presents 
a classification using conventional and topological order parameters. This classifi- 
cation offers a valuable perspective on a variety of conducting correlated electron 
systems, from the cuprate superconductors to the heavy fermion compounds. Con- 
nections are drawn, and distinctions made, between collinear/non-collinear magnetic 
order, bond order, neutral spin 1/2 excitations in insulators, electron Fermi surfaces 
which violate Luttinger's theorem, fractionalization of the electron, and the fraction- 
alization of bosonic collective modes. Two distinct categories of Z2 gauge theories 
are used to describe the interplay of these orders. Experimental implications for the 
cuprates are briefly noted, but these appear in more detail in a companion review 
article (S. Sachdev, cond-mat/0211005). 



1 Introduction 



The foundations of solid state physics reside on a few simple paradigms of 
electron behavior which have been successfully applied and extended in a wide 
variety of physical contexts. The paradigms include the independent electron 
theory of Bloch, its more sophisticated formulation in Landau's Fermi liquid 
theory, and the Bardeen-Cooper-Schrieffer (BCS) theory of electron pairing 
by an instability of the Fermi surface under attractive interactions between 
the electrons. 

In the past decade, it has become increasingly clear that these paradigms 
are not particularly useful in understanding correlated electron systems such 
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as the cuprate superconductors and the "heavy fermion" compounds. Many 
electrical and magnetic properties of these materials are rather far removed 
from those of the Fermi liquid, and remain poorly understood despite much 
theoretical work. The strong correlations between the electrons clearly makes 
the Fermi liquid an inappropriate starting point for a physical understanding 
of the many electron ground state. 

One approach to the strong correlation problem, advocated here, is to begin 
at the point where the breakdown of the Bloch theory is complete: in the 
Mott insulator. These are materials in which Bloch theory predicts metallic 
behavior due to the presence of partially filled bands. However, the strong 
Coulomb repulsion between the electrons leads to dramatically different insu- 
lating behavior, often with a very large activation energy towards conduction. 
A key ingredient in our discussion here will be a classification of such Mott 
insulators: we will characterize the ground states of Mott insulators with a 
variety of "order parameters" . Some of these orders will have a conventional 
association with a symmetry of the Hamiltonian which is broken in the ground 
state, but others are linked to a more subtle 'topological' order [1]. 

An implicit assumption in our approach is that these same order parameters, 
or their closely related extensions, can also be used fruitfully in a description 
of other correlated systems, which may be either metallic or superconducting. 
In the latter systems, it is sometimes the case that the true long-range order 
characterizes only a proximate Mott insulator, and the order parameter 'fluc- 
tuates' at intermediate scales. In such a situation, the ground state and its 
elementary excitations are adiabatically connected to the Bloch/BCS states, 
but the order of the proximate Mott insulator is nevertheless important in 
understanding experiments, and especially those that explore correlations at 
mesoscopic scales. The theory of quantum phase transitions offers a powerful 
method for controlled predictions for such experiments: identify a quantum 
critical point associated with the onset of long-range order, and use it to ex- 
pand back into the region with fluctuating order. The reader is referred to 
Ref. [2] for the application of such an approach to recent experiments in the 
cuprates. 

With this motivation, let us turn to the central problem of the classification 
of Mott insulators. Newcomers to the Mott insulator problem should consult 
Section III of Ref. [2] at this point for an elementary introduction to the 
microscopic physics. The most important degrees of freedom in Mott insulators 
are the quantum spins S a (r), a — x,y, z, which resides on the sites, r, of some 
lattice. The states on each site transform under the spin S representation 
of SU(2) (we are usually interested in S = 1/2), and the sites are coupled 
together with the Hamiltonian 



H= ]T J(r,r')S a (r)S a (r') + ..., 



(1) 



where J(r, r') are short-range exchange interactions (usually all positive, real- 
izing antiferromagnetic exchange), and the ellipses represent possible multiple 
spin couplings, all of which preserve full SU(2) spin rotation invariance. Here, 
and henceforth, there is an implied summation over repeated spin indices. 

At the outset, it useful to distinguish Mott insulators by whether their ground 
states break the SU(2) spin rotation invariance of the Hamiltonian or not. The 
paramagnetic states have (S a (r)) = and preserve SU(2) spin rotation invari- 
ance, while the magnetically ordered states we consider break spin rotation 
invariance in the following simple manner 

(S a (r)} = Re [$ Q e lK - r l ; (2) 



here K is the ordering wavevector, while & a is a three-component complex 
order parameter; this order parameter transforms as a vector under spin ro- 
tations, while $ a — > e~* Ka $ Q under translations by a Bravais lattice vector 
a. For simplicity we only consider systems with a single ordering wavevector, 
although the generalization to multiple wavevectors is not difficult. 

The quantum theory of the magnetically ordered states has been well estab- 
lished for a long time: one considers slowly varying quantum fluctuations of 
the field $ a in spacetime, and so obtains the Goldstone spin-wave excita- 
tions. Here, we wish to push the logic of this approach further: what happens 
when the quantum spin-wave fluctuations become strong enough to destroy 
the magnetic order in the ground state, and we reach a paramagnetic state 
with ($ a ) = ? We can think of this paramagnetic state as one with a 'fluctu- 
ating' $ Q order: does this tell us anything about the physical properties of the 
paramagnetic state ? We argue here that a surprising amount of information 
can gleaned from this seemingly naive approach, and it ultimately shows the 
way to a classification of both the magnetically ordered and the paramagnetic 
states. 



1.1 Landau theory 



We begin with the simple canonical procedure of considering the effective po- 
tential for $ Q fluctuations; our procedure here is general enough that it applies 
equally to both quantum and thermal fluctuations in insulators, metals, and 
superconductors. The structure of this effective potential is constrained by 
spin-rotation invariance and Bravais lattice translational symmetries, and the 
following low order terms are always allowed: 

u v 

V(Q a ) = S $ * a <S> a + 2*W«*/J + 2*aK*fi*P (3) 
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here s, u, v are phenomenological Landau parameters (u > 0, v > —u). If the 
value of K is commensurate with a reciprocal lattice vector, then additional 
low order terms can appear in the effective potential, but we defer considera- 
tions of such terms to later in this subsection. 

As is usual, we begin with a minimization of V over the values of the 3 com- 
plex numbers $ Q . For s > 0, we obtain the optimum value $ a = 0, which 
obviously corresponds to the paramagnet. For s < 0, we obtain two distinct 
classes of minima, which are not related to each other by any symmetry of the 
Hamiltonian. They are 



(A) Non-collinear spins, v > : $ a = n la + m 2a 

with rii^a real, ni a ri2a = and n\ a = n\ a = ——. 

Zi Li 

(B) Collinear spins, — u < v < : $ a = e td n a 

with n a real and n 2 = . (4) 

a u+v v ' 

From (2) we can see easily that in case (B) the average values of the spins at 
all r are either parallel or antiparallel to each other, while in (A) the average 
spins values map out a circular spiral. 

All the solutions in (A) or (B) represent physically distinct magnetically or- 
dered ground states, but the states within a category are degenerate and can 
be transformed to each other by symmetries of the Hamiltonian. It is useful to 
more carefully specify the manifold of degenerate magnetically ordered ground 
states. For (A), the ground state manifold is easier to decipher by a different 
parameterization of n lt2a which solves the constraints in (4): 

nia + in 2a = t ac z c (j* b z b , (5) 



where a, b extend over |, j, a a are the Pauli matrices, e a b is the antisymmetric 
tensor, and z a is a two-component complex field with \z^\ 2 + \ zi\ 2 = ^ — s/2u; 
note that z a transforms like a S — 1/2 spinor under spin rotations, while 

z a - e-^ 2 z a , (6) 



under translation by a Bravais lattice vector a. It is easy to check that (5) 
is in fact the most general solution of the constraints for case (A) in (4), 
but with a two-fold redundancy: z a and — z a correspond to the same non- 
collinear ground state. The two complex numbers £1X6 equivalent to four 
real numbers, and hence the manifold of ground states is S3/Z2, where S3 is 
the three-dimensional surface of a sphere in four dimensions. Turning to (B), 
n a maps out the surface of an ordinary sphere, S2, while the phase factor 
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e %e is U(l) = Si, a circle. However the factorization of $ a into n a and e %e is 
redundant because we can map n a — > — n a and — > + 7r without changing 
$ a ; hence the manifold of ground states for (II) is (S 2 x Si)/Z 2 . Summarizing, 
we have 



(A) Non-collinear spins, ground state manifold = S 3 /Z 2 ] 

(B) Collinear spins, ground state manifold = (S 2 x Si)/Z 2 . (7) 



This is a good point to mention additional restrictions that are placed on the 
ground state manifold at special commensurate values of K. These arise for K 
such that 2pK. equals a reciprocal lattice vector, where p is an integer. Then 
from the transformation of $q, under Bravais lattice translations we observe 
that the effective potential can contain the additional term 

Komm = -W (<5> a <5> a f + C.C. (8) 



For case (B) it is easy to see from (4) and (8) that V comm acts only on the 
angular field 9, and selects 2p values of 9 in the ground state. So the order pa- 
rameter manifold is now reduced to (£2 x Z 2p )/Z 2 : this change in the manifold 
has important consequences for p = 1, but is not of particular importance at 
larger values of p. For case (A), it follows from (4) and (8) that V comm has no 
influence on the ground state energy, and that the order parameter manifold 
remains S 3 /Z 2 . 

Having laid the ground work, and now are ready to extend this simple un- 
derstanding of magnetically ordered states to paramagnetic phases, where the 
order parameter $ a , constrained as in (4), is fluctuating. As we will see, global 
aspects of the order parameter manifold in (7) will play a crucial role. We will 
consider case (A) with non-collinear spins in Section 2, while case (B) with 
collinear spins will be discussed in Section 3. In both cases we will attempt 
to understand a variety of phases, including insulators, metals, and supercon- 
ductors. The evidence so far indicates that the cuprate superconductors are 
in category (B), and this is discussed at some length in Ref. [2]. Some of the 
more exotic phases appear in category (A), and we anticipate these will find 
realizations in heavy-fermion compounds, and particularly in those in which 
the magnetic moments reside on frustrated lattices such as the pyrochlore or 
the triangular. 
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2 Noncollinear spins 



An understanding of $ Q fluctuations in the paramagnetic phase requires that 
we proceed beyond the simple effective potential in (3), and consider spacetime 
dependent fluctuations using a suitable effective action. We assume that the 
non-collinearity of the spin correlations is imposed at some short scale, and 
so at longer scales we wish to impose the constraints under (A) in (4) at 
the outset in our effective action. This is most conveniently done using the 
parameterization in (5). We discretize spacetime on some regular lattice of 
sites j on which we have the spinors Zj a : note that this lattice may have little 
to do with the lattice of sites on which the underlying spins reside, and that we 
are working here on a coarse-grained scale. Any effective action on this coarse- 
grained lattice must be invariant under global spin rotations, and also under 
the global lattice transformation (6). However, most importantly, we note from 
(5) that it must also be invariant under the Z 2 gauge transformation 

z ja ~^ £ j z ja (9) 

where Sj = ±1 is a spacetime dependent field which generates the gauge 
transformation. This transformation is permitted because the local physics 
can only depend upon the order parameter which is invariant under (9). 

We begin in Section 2.1 by considering a very simple model of zj a fluctuations. 
This model is surely an oversimplification for the complex quantum systems of 
interest here, but it will at least allow us to obtain an initial understanding of 
possible phases and the global structure of the phase diagram. We will discuss 
applications to realistic physical systems in Section 2.2. 

2.1 Simple effective action 

Our model here omits long-range interactions and mobile excitations that can 
carry charge, which implies that it could apply directly only to insulators. We 
also neglect all Berry phases associated with the underlying spins. This will 
allow us to at least obtain a first understanding of the possible phases and 
the global structure of the possible phase diagrams. More careful considera- 
tions, described partly in Section 3.1 below, show that the Berry phases are 
crucial for realizing that 'bond order' is present in the 'confining' phase to be 
discussed shortly, but that they can be safely neglected in the other phases. 
The presentation below is based upon results contained in Refs. [3,4,5,6]. 

With the above caveats, we initially introduce the following simple partition 
function: 
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Z z = / Y[dz ja dz* ja 5 (\z ja \ 2 - l) exp (S z ) 
J j 

$z = -Jl Z ia Z ib Z jb Z ja -J2YI ( Z ia Z ib Z jaZjb + CC.) + . . . , (10) 

{ij) {ij) 

where we have rescaled the Zj a fields so that they obey the unit length con- 
straint on every site j, and (ij) represents nearest neighbors. The terms shown 
in (10) constitute the most general coupling between two sites consistent with 
the symmetries discussed above. Note that at this order the discrete symmetry 
(6) has effectively restricted us to terms which are also invariant under the 
U(l) symmetry Zj a — > e w Zj a with ip arbitrary: so the model S z has a global 
SU(2)xU(l) symmetry. 

While it is possible to proceed with (10), the full spectrum of possible phases 
is seen more easily by a representation which makes the Z 2 gauge invariance 
more explicit. For this, we must introduce a Z 2 gauge field = ±1, which 
resides on the links of the lattice, and which obeys the mapping 

O ij — ► EiGijEj (11) 

under the gauge transformation in (9). We can loosely view as arising from 
a Hubbard-Stratonovich decoupling of the quartic terms in (10). With in 
hand, we can now propose the alternative action as Z 2 gauge theory 



^noncoll ) 

J2Y Z ia Z ib Z jb Z ja- (12) 

{ij) 

The second term is the sum of the products of Uij over all elementary plaque- 
ttes of the lattice — it is the standard Maxwell term for a Z 2 gauge field. The 
first term proportional to J\ is easily seen to have a global 0(4) symmetry of 
rotations on the order parameter manifold S3. This symmetry is not present 
in the underlying spin model, and so we have added the last term propor- 
tional to J 2 which reduces the symmetry down to the required SU(2)xU(l). 
For K = 0, we can freely sum over the independently on each link, and 
the resulting action is easily seen to have a structure identical to S z in (10). 
For large K, the action (12) will allow us to easily access states which would 
have been harder to extract from (10). 

What is the physical significance of the flux in the Z 2 gauge field which is 
controlled by the coupling K in (12) ? It is a measure of the location of 
defects associated with the homotopy group tti(Ss/Z 2 ) = Z 2 , now often called 
'visons' [7]. Upon encircling a vison in a closed loop, the values of z a change 



^noncoii = / Y[dz ja dz* ja 5 {\z ja \ 2 - l) ^ exp (- 

Snoncoll = ~^ Z ia a 'ij Z ja + CC. - K ^ J| <7jj 

{ij) D {ij)ea 
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Fig. 1. A vison defect associated with tti(Ss/Z 2 ) = Z 2 . The left hand side represents 
a two-dimensional section of spacetime, while the right hand side is the space of 
points z a which reside on 53. Polar opposite points on 53 are represented by A and 
B. The Z 2 gauge field <7y = 1 along the loop, but = — 1 on the link directly 
connecting A and B. There is a Z 2 gauge flux at the center of the vison, and this 
flux allows identification of the vison defect even when the z a are fluctuating, and 
there is no magnetic order. 

smoothly along the loop, but the initial and final values of z a reside on polar 
opposite points on S3] all physical energies are independent of this sign change, 
and so there is no practical significance to this cut (see Fig 1). From (12) we 
see that the action will be minimized by = 1 along the loop, but we prefer 
Oij = — 1 across the cut between the initial and final points . Measuring the Z 2 
flux associated with this configuration of o^-, we conclude that the flux resides 
in a small region at the center of the vison. 

We can now sketch a phase diagram of iS non coii m the Ji,K plane in 2+1 
spacetime dimensions (related considerations apply to higher spacetime di- 
mensions). Many important features follow immediately by an analogy with 
a simpler problem considered by Lammert et al. [8,9] in the entirely different 
context of liquid crystals: they examined a problem with an order parame- 
ter belonging to the space S 2 /Z 2 (in contrast to the S 3 /Z 2 order parameter 
of interest here), which also permits an effective action in a Z 2 gauge the- 
ory very similar to <S n0 ncoii- Using their results, and those of the earlier work 
of Refs. [10,11], we obtain the phase diagram for 2+1 spacetime dimensions 
shown in Fig 2. 

Consider first the physics near the line K — 0. Here we can simply sum over 
the <7y, and work instead with the original action S z in (10). At large J\ we 
have a conventional non-collinear magnetically ordered phase with (z a ) 7^ 0, 
and hence 7^ 0. With decreasing Ji, there is a phase transition to a 

paramagetic phase. The key to understanding the nature of this paramagnetic 
phase, and of the quantum phase transition, is the observation that S z contains 
only gauge-invariant bilinears in Zj a on each site, and so we should be able to 
describe the physics in terms of the bilinear fields ni )2a in (5). Alternatively 
stated, the strong Z 2 gauge fluctuations (and the accompanying proliferation 
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Non-collinear magnetic order 



B 



Confining 
paramagnet, 
S=1 excitations 



Paramagnet with 
visons suppressed, 
S=1/2 excitations 











K 



Fig. 2. Schematic phase diagram of the model 5 nonco n in (12) as a function of J\ and 
K in 2+1 spacetime dimensions. We assume that J2 takes a small positive value. The 
model 5 non coii describes some of the physics of Mott insulators with non-collinear 
spin correations. However it omits Berry phase terms which impose the Heisenberg 
commutation relations of the spins in (1). These Berry phases play an important 
role in the small K region of the phase diagram: they induce bond order in the 
confining paramagnet, along with additional intermediate phases. The influence of 
the Berry phases is studied in Section 3.1 and in Ref. [12]. The suppression of visons 
in the deconfined state at large K is interpreted as a 'topological' order. 



of visons) has confined the z a quanta, and the elementary excitations of the 
paramagnetic state are the rti^ a quanta. Notice that these fields carry spin 
5 = 1, and so the paramagnetic phase is a confining state with only integer 
spin excitations. The quantum transition between the magnetically ordered 
state and this confining state, along the line AM in Fig 2 can be addressed by 
continuum models expressed directly in terms of the ni^a- This has been done 
by Kawamura [13], Pelissetto et al. [14,15], and Calabrese et al. [16] in the 
context of classical phase transitions in stacked triangular lattices, and they 
obtained a continuous phase transition in the "chiral" universality class in 
2+1 spacetime dimensions. The global symmetry at this "chiral" fixed point 
remains SU(2)xU(l). 

Let us turn now to the physics near the line J\ = in Fig 2. Here, the 
z a fluctuate strongly and can be integrated out, leading to a pure Z 2 gauge 
theory in the (Ty. In 2+1 spacetime dimensions, this gauge theory is well 
known [10,11] to have a confinement- deconfinement transition at some critical 
K near the point C: the large K region is the deconfined state where vison 
fluxes are suppressed. Consequently, in this deconfined state, we can use a 
simple physical picture in which we choose a gauge with all = 1: the 
z a quanta in (12) are then free to propagate through the system. The large 
K, small J\ region is therefore a paramagnetic state with neutral 5 = 1/2 
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'spinon' excitations. The suppression of visons in this state is interpreted as 
the presence of 'topological order'; this is similar to the topological order below 
the Kosterlitz-Thouless transition in the two-dimensional classical XY model, 
where point vortices are suppressed [1]. 

Finally, to complete our picture of the phase diagram, we look at the region of 
large K. Here, visons are suppressed, and so we can freeze = 1, and we are 
left with a simple model of interacting z a quanta, with a global SU(2)xU(l) 
symmetry. Such a theory has been studied in its continuum limit, and yields a 
description of the magnetic ordering transition along the line BM: the critical 
point has a large 0(4) symmetry in 2+1 spacetime dimensions, and the critical 
exponents are those of the 4-component tp A field theory [17,5]. 

2.2 Physical Applications 

The results of Section 2.1 for large K apply directly to Mott insulators with 
non-collinear spin correlations. The mapping of the magnetically ordered phase 
is evident, while we identify the paramagnet with visons suppressed with the 
resonating valence bond (RVB) state of Refs [18,19,20,21,22,23]. The latter 
identification follows from the full preservation of symmetries of the Hamil- 
tonian, the presence of neutral S — 1/2 spinon excitations, and from the 
topological order implied by the suppression of visons [3,22,7]. Experimen- 
tally, we note that the evidence for a RVB state in Cs 2 CuCi4 [24] is in a 
system with non-collinear spin correlations, consistent with our theoretical 
picture. At small K, and smaller J 1; the proliferation of visons requires more 
care, as the quantum mechanics of the underlying spins contributes a Berry 
phase to each vison [25,7]. These phases induce bond order in the confining 
paramagnet, and can also lead to additional intermediate states [7,12]: we will 
not discuss this further here, but will illustrate the influence of Berry phases 
in the simpler context of collinear spins in Section 3.1. 

Let us finally move away from the Mott insulator, and introduce mobile charge 
carriers into the RVB state just discussed, by doping in holes. A scenario 
which has been intensively discussed in the literature [20,21,7] is that hole 
fractionalizes into independent quasiparticle excitations: a neutral S — 1/2 
'spinon' and a spinless, charge e holon. The spinon here is, of course, that 
just discussed above here in the proximate Mott insulator. We can express 
this fractionlization of the injected hole by the following schematic operator 
relation for the electron annihilation operator c a : 

Ca = f (AiZ + \2tab4) » ( 13 ) 

where z a is the bosonic neutral spinon operator from Section 2.1, p is a 
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fermionic operator that creates a spinless hole, and Ai ; 2 are constants that de- 
pend upon the details of the microscopic physics. It is possible for the fermionic 
statistics to pass from the hole to the spinon by binding between quasiparticles 
and visons, as has been discussed in some detail in Refs. [26,27,28]. Then the 
relationship (13) would be replaced by 



where f a is a neutral, S = 1/2, fermionic spinon, b is charge e bosonic holon, 
and A3 5 4 are constants. In this scenario, the Bose condensation of the holon 
leads to a BCS superconductor with vestiges of the topological order of the 
RVB Mott insulator: experimental probes of this 'fluctuating' topological order 
have been proposed [29,30,31,2], but no positive signal has been observed so 
far in the cuprates [32,33]. 



2. 2. 1 Fractionalized Fermi liquids 

It has been argued recently [34,35] that a different possibility is more likely 
for the doped RVB Mott insulator in d > 2 spatial dimensions. The doped 
electrons (or holes), instead of fractionalizing into spinless charged particles 
and neutral S — 1/2 spinons, retain their integrity in the ground state, and 
their spin and charge remain bound to each other. At the same time, the 
neutral S — 1/2 spinons of the RVB Mott insulator survive in the doped 
system, and are only renormalized slightly by the mobile carriers. So, alter- 
natively stated, the added electrons form a Fermi-liquid-like state which is 
approximately decoupled from the spinons. The resulting metallic state has a 
Fermi surface with S = 1/2, charge — e quasiparticle excitations, along with a 
separate set of neutral spinon excitations which continue from the RVB Mott 
insulator. Moreover, the volume enclosed by the Fermi surface is 'small' i.e. it 
is determined solely by the density of doped electrons, and does not include 
the spins of the Mott insulator. In this situation, the volume of the Fermi sur- 
face becomes a direct experimental signal of the topological order in the RVB 
state. Heavy fermion compounds on frustrated lattices with weak or absent 
magnetic order are likely candidates for realizing this state. 

The fractionalized Fermi liquid defined above violates the standard Luttinger 
theorem, and it is useful to express this violation in its most general terms. 
Consider a correlated electron system on a periodic lattice in d spatial dimen- 
sions, whose ground state preserves time-reversal and spin rotation invariance, 
and has unit cell volume vq. Let ut be the total density of electrons per vol- 
ume v ; tit includes all the electrons in the system, including e.g. the core Is 




(14) 
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electrons. Luttinger's theorem states that in a conventional Fermi liquid state 
2 x x (Volume enclosed by Fermi surface) = n^mod 2). (15) 

{27l) d 



The leading factor of 2 on the left hand side comes from spin degeneracy, while 
modulus 2 on the right hand side allows fully filled bands to not contribute 
to the Fermi surface volume. The fractionalized Fermi liquid being discussed 
here violates (15), but instead obeys a modified relation 

2 x V ° x (Volume enclosed by Fermi surface) = (n T — l)(mod 2). (16) 
(2ir) d 



So exactly one electron per unit cell has dropped out from the Fermi volume. 
This can only happen in a topologically ordered state which possesses neu- 
tral S — 1/2 spinon excitations, in addition to the electron- like quasiparticle 
excitations on the small Fermi surface. 



3 Collinear spins 



We turn here to the second broad category of correlated electron systems 
introduced in Section 1: the collinear spin case (B) in (4). We will find that the 
paramagnetic states with fluctuating collinear spin order are entirely different 
from those present for non-collinear spins in Section 2.2. The order parameter 
manifold, from (7) is now (S 2 x Si)/Z 2 , and the appearance of a Z 2 quotient 
means that the physics can again be described in a generalized phase diagram 
of a Z 2 gauge theory [36]. However, the Z 2 flux now identifies a new type 
of defect which we will discuss in Section 3.2. The Z 2 gauge theory also has 
a deconfined phase, but this Z 2 fractionalization does not lead to neutral 
S — 1/2 excitations; instead, as we will see in Section 3.2 it is the spin and 
charge collective modes which 'fractionalize' apart from each other. 

A separate crucial property of Mott insulators with collinear spin correlations 
in spatial dimension d — 2 is the ubiquity of confining states with only integer 
spin excitations [37]. Moreover, except for certain special values of the spin S 
per unit cell, these confining states break lattice space group symmetries by the 
appearance of spontaneous 'bond order' in the ground state. This bond order 
also often survives in proximate conducting states obtained by doping the Mott 
insulator [4,38]. Let us define bond order more precisely here: most generally, 
bond order implies a modulation in observables invariant under spin rotation 
and time reversal which break a space group symmetry of the Hamiltonian. 
The simplest such observable in states proximate to Mott insulators is simply 
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the spin exchange energy: so we can introduce the bond order parameter 
Q a (r) = (Sa(r)S a (r + a)>, (17) 



where a is usually a vector connecting near neighbors. Note that for a = 0, 
Qo( r ) is a measure of the charge density on site r, and so (if permitted by the 
symmetry of the state) there is usually also a modulation of the site charge 
density in a bond-ordered state. However, we expect the long-range Coulomb 
interactions to suppress modulations at a = 0, while those with a^O should 
be significantly larger. 

We will begin our discussion in Section 3.1.1 by a detailed discussion of the 
simplest, and most common, Mott insulator with collinear spin correlations: 
that on a rf-dimensional cubic lattice with ordering wavevector K = (jr, ir, . . .). 
This order is commensurate with p — 1, in the notation of (8). We will explic- 
itly evaluate the quantum spin Berry phases for this case, show their intimate 
connection to bond order. Section 3.1.2 will also include extensions of the 
study of bond order to doped systems with mobile carriers. 

States at other wavevectors and in metals and superconductors, and the Z2 
gauge theory of the fractionalization of their excitations will be discussed in 
Section 3.2. 



3.1 Berry phases and bond order 



3.1.1 Mott insulators 

Consider the insulating antiferro magnetic (1) on a <i-dimensional cubic lattice 
with predominant nearest-neighbor exchange interactions. This should have 
collinear spin correlations with K = (n, ir, ...), which allows the term with 
p — 1 in (8). Inserting <3> a = e t6 n a (obtained from (4)) into (8) we observe 
that we can always choose the origin of co-ordinates so that the values 6 = 0, 71 
are preferred. These do not lead to new values of $ Q , and so in this case we 
have simply $ a = n a , a real three-component vector. 

Now express the coherent state path integral of (1) using the values of the 
field n a on a d + 1 dimensional hypercubic lattice discretization of spacetime. 
The derivation of this path integral is reviewed in Chapter 13 of Ref. [39], and 
leads to the following partition function 

Zn = Jl[ dn i^ (rfa - 1) exp ( n ia n ja -iSj^ Vj^jr J , (18) 

3 V 3 fo) 3 J 



13 



where i,j are sites of the d + 1 dimensional hypercubic lattice (the symbol 
% = V - 1 i n the prefactor of the second term, and the context should prevent 
confusion on its meaning), and we have rescaled all the rij a to make them unit 
length. The first term in the action above is the analog of the terms in (10), 
and imposes a cost in the action for deviations from the perfectly ordered 
state; we expect a magnetically ordered state (n a ) ^ for small values of the 
coupling g, and a paramagnetic state with (n a ) = for large g. The second 
term in (18) is the all important Berry phase: here rjj = e lK ' rj = ±1 is a fixed 
field identifying the spatial sublattice of the site j (note that rjj is independent 
of the imaginary time co-ordinate, r). Finally Aj^, with \i = r, x, y, . . . taking 
d + 1 dimensional possible values, is defined by 

Aju = oriented area enclosed by the spherical triangle with vertices 

n j a , n 

j+fi,a, and any fixed reference N 0a , (19) 

where j + p, is the nearest site to j in the \x direction. It is customary to choose 
No a = (0,0, 1), the north pole, but it is not difficult to see that the partition 
function Z n is independent of the value of N 0a . Indeed, elementary geometric 
considerations of spherical triangles show [12] that a choice of a different A^ a 
leads to A'^ which is related to A^ by a 'gauge transformation' 

A' jlx = A m - A^j (20) 



where A M is the discrete lattice derivative in the \x direction, and (f>j is the area 
of the spherical triangle formed by nj a , N 0a and Nq . It should also be noted 
here that the area of a spherical triangle is uncertain modulo 4n, and so is the 
value of Aju, but the partition function Z n is insensitive to this uncertainty 
because e* 47r5 = I. 

For small g in Z n , fluctuations of n a are suppressed, and in dimensions d > 1, 
we are in the conventional Neel ordered ground state with (n a ) ^ 0. As an 
aside, we note the small g behavior for d — 1. In d — 1 we can continue to 
assume that n a varies smoothly from site-to-site for small g, and hence take 
the naive continuum limit of (18). This yields the 0(3) non-linear sigma model 
in 1+1 dimensions, along with a topological 6*-term with co-efficient 9 = 2irS, 
as reviewed in Ref. [39]. For integer S this exhibits the Haldane gap state, 
while for half-odd-integer S a gapless critical state as in the Bethe's solution 
of the nearest neighbor antiferromagnetic chain can appear. 

Our interest here is primarily in the larger g regime of Z n , where there are 
significant fluctuations of n a , and we are in a paramagnet with (n a ) = 0. 
There are large fluctuations here in the value of the Aj^ also, and so the 
Berry phase term requires careful evaluation. An explicit evaluation of this 
term is essentially impossible, but considerable progress has been made in the 
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'easy-plane' case, where spin-anisotropies reduce fluctuations to the equator 
in spin space with n,j a = (cos$j, sin9j, 0) [12]. Here we follow a simple 'hand- 
waving' argument [40] whose results are known to be consistent with all cases 
in which more sophisticated methods are possible; the paragraph below is 
partly reproduced from the review in Ref. [40], to which we will also refer the 
reader for additional details. 

The main idea for the larger g regime is to change variables from the order 
parameter rij a to the field which naturally represents the Berry phase. 
Formally, this can be done by introducing new 'dummy' variables A, M and 
rewriting (18) by introducing delta-function factors which integrate to unity 
on each link; this leads to 

Z n = f n^exp l-t2Sj2v 3 A JT ) I \{dn 3 J(n%-l)5{A m /2- A m ) 

Jl> V 3 / 3 

x exp I -^^Z n ^ n i 
= J n^>exp ( -S A {A^) - i2Sj2vjA jr ] • (21) 

J> \ 3 J 

In the first expression, if the integral over the A, M is performed first, we trivially 
return to (18); however in the second expression we perform the integral over 
the rij a variables first at the cost of introducing an unknown effective action 
Sa for the A 3il . In principle, evaluation of Sa may be performed order-by- 
order in a "high temperature" expansion in 1/g: we match correlators of the 
Aj^ flux with those of the Aj^ flux evaluated in the integral over the rij a 
with positive weights determined only by the 1/g term in (18). Rather than 
undertaking this laborious calculation, we can guess essential features of the 
effective action S A from some general constraints. First, correlations in the rij a 
decay exponentially rapidly for large g (with a correlation length ~ 1/ ln(g)), 
and so Sa should be local. Second, (20) implies that S A should be invariant 
under the gauge transformation 

A' m = A m - A^./2. (22) 



Finally, uncertainity of Aj^ modulo 4n implies that Sa should also be invariant 
under 

A m - A m + 2tt. (23) 



The simplest local action which is invariant under (22) and (23) is that of 
compact U(l) quantum 'electrodynamics 7 and so we propose that for larger 
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values of g, Z n Za with 

Z A = [ [] dA m exp f 1 E cos ( A » A i" ~ A » A in) ~ l2S E 5 ( 24 ) 

comparison with the large g expansion shows that the coupling e 2 ~ g 2 . This 
is one of the central results of this subsection [37,41,42,12]: the paramagnetic 
state of a Mott insulator with two-sublattice collinear spin correlations in d 
spatial dimensions is described by a compact U(l) gauge theory in d+1 space- 
time dimensions 1 , accompanied by a Berry phase term as specified in (24). 
In the gauge theory/particle physics language, the Berry phase corresponds 
to static 'matter' fields of charge ±2S residing on the two spatial sublattices 
of the c?-dimensional lattice. 

We are now faced with the technical problem of evaluating the partition func- 
tion of compact QED in d + 1 spacetime dimensions in the presence of static 
background matter fields, as written in (24). There is a good understanding 
of the phases of this theory in d — 1, 2, and the reader is referred to another 
review [40] for further details — here we will summarize the main conclusions. 
The structure in d = 3 is not yet fully understood, and this remains an impor- 
tant avenue for future research. All of these previous analyses rely on a duality 
mapping of (24); the mapping proceeds in a canonical manner, beginning with 
a rewriting of the cosine term in (24) in the Villain periodic Gaussian form, 
followed by the evaluation of the integral over the A, M , which results in a par- 
tition function of 'dual' fields. The advantage of this duality mapping is that 
the Berry phases are exactly accounted for at the outset, and the partition 
function in the dual representation has only positive weights — this makes it 
amenable to standard analyses of statistical field theory. A discussion of the 
results in various spatial dimensions, d, follows. 

In d = 1, the duality mapping of the Villain form of Za leads to a simple 
partition function that can be evaluated exactly [40,12]. There is a gap to 
all excitations, and for half-odd-integer S a two-fold degenerate ground state 
and a broken translational symmetry associated with the appearance of bond 
order, as illustrated in Fig 3. 

In d = 2, the duality mapping of (24) leads to interface models (also called 

1 Some readers may find the following alternative interpretation of the compact 
U(l) gauge theory useful. Express rij a = Wj a cr ab w jb where Wj a is a two-component 
complex spinor. This parameterization has a U(l) gauge redundancy corresponding 
to Wj a — > e^Wja (contrast this with the gauge redundancy associated with (5)). 
Expressed in terms of the Wj a , the theory is a lattice discretization of the CP 1 
model [43,44,12], with Aj^ its compact U(l) gauge field. Integrating out the Wj a in 
the paramagnetic phase leads to (24). 
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Fig. 3. Schematic representation of the ground state of for half-odd- integer S 
in d = 1. The state is identified with the different values of (5a( r ) in (17) with 
a a nearest neighbor vector: the links with ellipses have values of Q a (r) distinct 
from those without. The reader can also view this picture as caricature of the 
wavefunction — a simple trial state for S = 1/2 is a product of singlet valence bonds 
between sites linked by an ellipse 

height or solid-on-solid models) in 2+1 dimensions [37,42,12] (remarkably, the 
same interface models are also obtained by a duality mapping [41,45] on quan- 
tum dimer models [46] of the paramagnetic state). In statistical mechanics, 
these interface models would describe the growth of a three-dimensional in- 
terface of a four- dimensional crystal. The Berry phases in (24) lead to offsets 
in the allowed configurations of the interface model. A fundamental property 
of interface models in 2+1 dimensions is that their configurations are smooth 
for all couplings. In the smooth state, the symmetry of uniform translations 
of the interface by a constant is broken, and the interface has some fixed av- 
erage height. When combined with the offsets produced by the Berry phases, 
the smooth interface is seen to imply a broken lattice space group symmetry 
because of the appearance of bond order (except for certain special values of 
S which depend upon the lattice structure — for the square lattice, there is 
no bond order for even integer S). The fluctuations of the smooth interface 
are gapped, and these correspond to a gapped singlet excitation mode of the 
quantum antiferromagnet. In the QED langauge, the gauge-field excitations 
are gapped, and the compact U(l) gauge theory is in a confining state. The 
spinful excitations of the antiferromagnet consist of a gapped S = 1 mode 
associated with the fluctuations of n a — there are no neutral excitations which 
carry S = 1/2. 

A full analysis of the situation in d = 3 has not been carried out. However, it 
is clear that a new state can appear in the phase diagram-the compact QED 
theory has a gapless Coulomb phase with a 'photon' excitation [47,48,49,50]. In 
the quantum antiferromagnet, this photon corresponds to a collective spinless 
excitation with a gapless linear dispersion. It is likely that this Coulomb phase 
is insensitive to the Berry phases, has no bond order, and allows neutral S = 
1/2 spinon excitations. The analog of the bond-ordered states of d = 2 should 
also be present. 

Large scale numerical studies of S = 1/2 quantum antiferromagnets on the 
square lattice which allow for the destruction of magnetic order in the ground 
state have only just become possible. Two studies have appeared recently 
[51,52], and both present unambiguous evidence of bond order as shown in 
Fig 4. 
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(a) (b) 

Fig. 4. Sketch of the two simplest possible states with bond order for S = 1/2 on 
the square lattice: (a) the columnar spin-Peierls states, and (b) plaquette state. 
The different values of Q a (r), defined in (17), on the links are encoded by the 
different line styles. Both states are 4-fold degenerate; an 8-fold degenerate state, 
with superposition of the above orders, also appears as a possible ground state of 
the generalized interface model. 

3.1.2 Mobile charge carriers 

The fate of the d = 2 bond-ordered paramagnet upon doping with mobile 
charge carriers has been investigated in a number of studies [4,38,53,54,55]. 
As stated earlier, the bond order survives for a finite range of doping concen- 
trations. Moreover, the ground state also acquires <i-wave-like superconduct- 
ing order which co-exists with bond order, a state first predicted in Ref. [4]. 
Loosely speaking, this superconductivity can be understood as a consequence 
of the mobility of the singlet valence bonds in Fig 3: the mobile pairs of elec- 
trons behave as Cooper pairs, and their Bose-Einstein condensation leads to 
superconductivity. The configuration and period of the bond order can also 
evolve with increasing density of carriers [38,53,55], especially if the parame- 
ters are such that the physics of frustrated phase separation is important [56]. 
The reader is referred to Ref. [2] for a recent review of this work, along with 
implication for a number of recent experiments on the cuprate superconduc- 
tors. These experiments include a possible observation in STM of a state with 
co-existing bond order and <i-wave superconductivity [57,58,59,55,60]. 



3.2 Collective mode fractionalization 

We now consider magnetically ordered and paramagnetic states with collinear 
spin correlations at wavevectors other than K = (n,7r, ...). For these cases 
the field e td , in the decomposition $ Q = e l0 n a in (4), cannot be disposed 
at the outset, and needs to be retained as an additional degree of freedom. 
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The physical significance of this additional field is readily apparent from a 
comparison of the definitions (2) and (17) of the magnetic order parameters; a 
factorization of the expectation values in (17) implies that there is contribution 
[61,62]: 



Qa(r) = Re 



e 2i0 e iK-(2r+a) 



+ 



(25) 



where we have assumed, as in Section 3.1.1, that n a has been rescaled so that 
r? a — 1. So a collinear magnetically ordered state at wavevector K automat- 
ically has co-existing bond order modulations at wavevector 2K, and e 2ie is 
the bond order parameter. 

We can now combine this observation with the results of Section 3.1.1 to 
obtain a simple, and quite general, theory of the loss of magnetic order for the 
values of K under consideration here. As bond order co-exists with magnetic 
order, it is natural consider the possibility that the loss of magnetic order 
occurs in a 'background' of bond order, and the latter is present on both sides 
of the transition. In other words, we consider a transition from a magnetically 
ordered state with (3> Q ) 7^ (and so necessarily (e 2ie ) 7^ 0), to a paramagnetic 
state with bond order which has (<& a ) = but still (e 2ld ) 7^ 0. We can assume 
that the background bond order is quenched, and hence need only consider a 
theory of n a fluctuations in this environment. 

For insulating systems, this theory will have essentially the same structure 
as that in Section 3.1.1, with the difference that the coupling constants will 
now acquire a modulation induced by the background bond order. The Berry 
phases will again attempt to induce bond order in the regime with strong 
fluctuations of n a , but this may be a redundant effort-the bond order is already 
present in the underlying Hamiltonian for the spins. As long as the underlying 
bond order has an even number of spins per unit cell, the upshot is that we 
can safely neglect the Berry phases, and the theory is simply that of a n a field 
undergoing a 'quantum disordering' transition in a partition function with 
positive weights and short-range interactions: the critical theory for this is the 
0(3) ip 4 field theory in 2+1 spacetime dimensions. 

For systems with mobile charge carriers, we have to consider possible long- 
range couplings induced by other gapless excitations. These effects are quite 
strong and relevant in metals, but are relatively unimportant in the supercon- 
ducting states considered in Section 3.1.2: superconductors have low energy 
fermionic excitations at only select points in the Brillouin zone, and these 
usually couple quite weakly to the critical spin fluctuations [53]. 

Having disposed of this simple possibility for the loss of collinear magnetic 
order, let us now explore more fully the possible phases that are allowed in a 
general interplay of collinear magnetic and bond order parameters. The list of 
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fractionalization transition is between phases with and without suppression of 
these defects. Both phases have neither magnetic or bond order. 

possibilities is very large, but interplay with experiments should help narrow 
the range of possibilities [36,63,64]. As in Section 2.1 a global view of the 
phase diagram is obtained most easily by constructing an effective action for 
a Z 2 gauge theory. The gauge theory must now be formulated for an order 
parameter on the manifold (S 2 x Si)/Z 2 , rather than S 3 /Z 2 , but the strategy 
is quite similar to that of Section 2.1. We take the physical field <& a , and split 
it apart into its n a and e %e constituents; the overall sign for the constituents 
involves a Z 2 gauge choice which can be made differently at distinct points in 
spacetime. We can compensate for this ambiguity by introducing a Z 2 gauge 
field, <7jj (note that the physical interpretation of this is entirely different 
from that in Section 2.1), and so obtain the effective partition function [36]: 




For commensurate values of K there will be an additional on-site anisotropy 
field arising from V comm in (8) ~ — w J2j cos(2p9j). 

What is the interpretation of a Z 2 gauge flux now ? It locates the position 
of defects often called "stripe dislocations", sketched in Fig 5. To see this, 
consider, as in Section 2.1, the defects associated with tti((S 2 x Si)/Z 2 ) [36]. 
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These include ordinary vortices in the field 6, under which 6 winds by an 
integer multiple of 2tt upon encircling the vortex. However, an additional class 
of half-vortices are also allowed for which the Z 2 quotient plays a crucial role: 
in these, 9 winds by an half-odd-integer multiple of 2rr. At the same time, there 
is a corresponding trajectory for n a which connects polar opposite points on 
S 2 ; this ensures that the physical field $ a is single-valued everywhere. In the 
context of (26), we observe that, as in Section 2.1 and Fig 1, we will have 
o"jj = — 1 only between the initial and final points of the loop, to eliminate the 
branch cuts in e td and n a . Hence each such half vortex carries a Z 2 gauge flux. 

The phase diagram of S co \\ in the three-dimensional ^4, K space is quite 
complex, but many of its features can be understood by analyses of various 
two-dimensional sections, most of which bear some formal similarity to the 
model analyzed in Section 2.1. We do not want to run through the plethora of 
possibilities here, and will merely focus on a particular phase transition [63,36] 
which has some intriguing experimental implications. 

Consider the region of small J 3i4 where (§ a ) = and (e 2td ) = 0, and so there 
is no magnetic or bond order. Let us first take a small value of K — in the 
cuprates, we would imagine this is in the overdoped region; then the field 
will fluctuate strongly, and the stripe dislocations will proliferate. Such a state 
has no vestiges of either magnetic or bond order, and is therefore an ordinary 
superconductor or Fermi liquid. The strong o~ij fluctuations also 'bind' the 
n a and e l6 fields to each other in the <& a field, and the latter constitutes a 
collective elementary excitation which carries both spin and bond correlations; 
indeed, $ a is simply an ordinary collective particle-hole excitation which can 
be described by a time-dependent Hartree-Fock theory; in a superconductor, 
this excitation could also avoid damping from the particle-hole continuum, and 
so become a sharp quasiparticle. This collective mode has 6 real components, 
corresponding to the 3 complex numbers required to define $ a [65]. 

Now increase the value of K; in the cuprates, this would correspond to moving 
towards the underdoped region. This will suppress the stripe dislocations until 
eventually there is a Z 2 gauge theory transition to a deconfined phase; this 
transition is the analog of the transition across the line MC in Fig 2. The 
deconfined phase has dislocations suppressed, and unbound quanta of n a and 
e %e . However it is still true that = and (e 2td ) = 0, and so there is no 

magnetic or bond order. So this is a fractionalized phase in which quanta of 
$ a have separated into 3 real quanta of n a and 2 real quanta of e l6 i.e. the 
collective mode has fractionalized. Note that unlike the fractionalized phase of 
Section 2.1, there are no neutral S — 1/2 excitations, and the electron always 
retains its bare quantum numbers. 

Why is this Z 2 gauge fractionalization transition interesting ? First, unlike 
that in Section 2, it is based on correlations (collinear magnetic and bond) 
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that are actually observed in the lightly doped cuprates. While long-range 
order with 7^ and (e 2t6 ) 7^ is observed at low doping in some of 

the cuprates, it is quite clear that all such conventional order has essentially 
disappeared by the time we reach optimal doping. However, remnants of such 
order may still be present in that the stripe dislocations of Fig 5 have not yet 
proliferated. This then opens up the possibility [36,63] of a quantum critical 
point near optimal doping [66,67], associated with the proliferation of stripe 
dislocation with increasing doping, which is in the universality class of the 2+1 
dimensional Z 2 gauge theory. As an added bonus, it has been recently argued 
[68] that this gauge theory remains strongly coupled even in the presence 
of Fermi surface damping, and so is an attractive candidate for anomalous 
behavior at higher temperatures. 



4 Conclusions 

This article has given a broad overview of phases of strongly interacting elec- 
trons, going beyond those found in the conventional Bloch-BCS theory of solid 
state physics. A plethora of exotic and conventional states have been proposed 
in the literature, and we attempted to bring some 'order' into the subject, by 
presenting a classification based on the physics of Mott insulators, and by 
making connections between the correlations in different states. 

Some of the states considered here have conventional order parameters, and 
could, in principle, also be obtained from a weak-coupling instability of the 
Hartree-Fock/BCS theory of the Fermi liquid or superconductor. These include 
the states with collinear or non-collinear magnetic order, and with bond order. 
We have chosen here, instead, to motivate them from the states of the Mott 
insulator. Our approach properly accounts for the strongest energy scales at 
the outset (the local repulsion between nearby electrons), and so has a better 
handle on the energy scales of various excitations. The Hartree-Fock/BCS 
theory could obtain states with the same overall symmetry, but would give very 
different estimates of the relative stability and energies of various excitations. 
The selection of the wavevector, K, would also be based upon the Fermi 
surface geometry. In contrast, the Mott-insulator-based approach used here, 
selects wavevectors from very different criteria: e.g. in the bond ordered states, 
it is the spin Berry phases and the resonance and alignment of singlet valence 
bonds which selects the state, and leads to predictions for the values of K 
[4,38] which have little to do with the Fermi surface geormetry. These latter 
values appear to be consistent with observations [69]. 

Another value in starting with the Mott insulator is that it also allows a 
systematic classification of states with more exotic order parameters, which 
cannot be expressed in terms of local correlation functions. They appeared 
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here as 'quantum disordered' states of different types of magnetic order; these 
'disordered' states however possessed a topological order associated with the 
suppression of certain defects which were classified by the global topology of 
the magnetically ordered state. We also found a description of the dynamics 
of this topological order in two distinct Z2 gauge theories, in systems with 
non-collinear and collinear spin correlations respectively. 

With some understanding of global features of the phase diagram at hand, 
we can proceed with a study of possible quantum critical points, and their 
influence on crossovers at finite temperatures [39]. For the cuprates, these 
issues have been reviewed in Ref. [2]; the 'stripe fractionalization' transition 
of Section 3.2 is only poorly understood, and much additional work is required 
to understand its observable consequences. The most intensive investigation 
of quantum criticality has so far occurred in the heavy fermion compounds, 
and most observations do not agree with the weak-coupling spin density wave 
picture. Here, we think that the small Fermi surface state of Section 2.2.1 
offers the prospect of leading to many interesting related states, and of non- 
trivial quantum critical points between them — this is an exciting avenue for 
further research. 
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